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Fibonacci system with both diagonal and off-diagonal dis- 
order is mapped to a model with dimer type 'defects'. The 
model exhibit resonance transition corresponding to zero re- 
flectance condition resulting in extended states in the system. 
These Bloch-type states reside on the line of inversion sym- 
metry of the energy spectrum and include states at the band 
crossings. The resonant states at the band crossings are re- 
lated to the harmonics of the sine wave and have transmission 
coefficient equal to unity in quasiperiodic limit. This is in 
contrast to other resonant states where the transmission co- 
efBcient oscillates. An exact renormalization scheme confirms 
the fact that all resonant states are Bloch type waves as they 
are described by trivial attractors of the renormalization fiow. 

PACS numbers: 71.23-k-h61.44-n 



Fibonacci systems have emerged as an interesting class 
of quasiperiodic systems exhibiting self-similar fractal 
wave functions as well as Bloch-type extended states. 

In earlier studies interest in these systems was 
due to the fact that they exhibit cantor spectrum and 
power-law decaying wave function and hence are in be- 
tween periodic and random systems. Recent interest in 
the Fibonacci systems is due to the possibiliy of extended 
states as was first discussed briefly by Wu et al It 
was argued that the Fibonacci systems with both diag- 
onal and off-diagonal quasiperiodic order are similar to 
the the binary alloy model which is dual to the dimer 
model. [|j Therefore Fibonacci system should exhibit the 
resonance predicted for the binary alloy model. Recently, 
Macia et al Q discussed the existence of extended states 
in these systems by imposing the condition of commuta- 
tivity of two non-equivalent blocks of transfer matrices. 
They obtained an explicit expression for the transmission 
coeflticient for a flnite size Fibonacci system and charac- 
terized these states as extended critical due to the fact 
that in the quasiperiodic limit the states have transmis- 
sion coefficient (T) which oscillates as the size of the sys- 
tem is varied. They argued that the quasiperiodic order 
play a key role in determing the extended states which is 
in constrast to the analysis of Wu et al where the exis- 
tence of extended states is due to a resonance condition 
unaffected by the long range quasiperiodic order in the 
system. 

In this paper, the Fibonacci lattices with both diago- 
nal and off-diagonal disorder are mapped to the lattices 
with dimer-type defects. Therefore, our results put the 
Fibonacci systems on the same footing as the lattices 
with dimer type correlations |q] and hence, the origin of 



extended states in this system is traced to the resonance 
transition corresponding to a destructive interference of 
the reflected wave between the two neighboring sites of 
the dimer defect. Using this mapping, we derive the reso- 
nance condition and establish the fact that the extended 
states as discussed by Wu et al are identical to the propa- 
gating states discussed by Macia et al. The fact that 
the extended states of the Fibonacci system are in fact 
Bloch waves on a decimated lattice is further conflrmed 
by an exact renormalization group (RG) scheme. 

Novel result of our studies is the relationship between 
the resonant states and the states where the energy bands 
cross. We show that the resonant states at the band 
crossings are fully transmitting states (T = 1) in the 
quasiperiodic limit and are described by the wave func- 
tions that are related to the harmonics of the sine wave 
with fundamental Bloch number equal to the golden 
mean. Therefore, unlike the rest of states with oscil- 
lating T, these states have a special property namely for 
Fibonacci sizes the perfect transmission limit is well de- 
flned. 

We consider a tight binding model (TBM) with both 
diagonal and off-diagonal disorder. 
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The Fibonacci system that we study is obtained by gen- 
erating a Fibonacci sequence from two symbols s and h 
by the substitution s b and b ^ bs. The corresponding 
onsite potential takes two values and £(,• The resulting 
lattice can be viewed as made up of bbs and sbs types of 
sub blocks. The off-diagonal couplings are of two types 
corresponding to two possible nearest-neighbors and are 
denoted as t^s and i^f,. We choose if,s = 1 and tf,b = 7 as 
shown in flgure 1. 

In order to map the Fibonacci system to a dimer 
model, we decimate all the s sites of the lattice with 
onsite potential eg. The remaining 6 type sites renor- 
malize in two different ways depending upon whether it 
belongs to bbs block or sbs block. As shown in flgure 1, 
the resulting decimated lattice can be viewed as a pure 
r type lattice with d type defects that always appear in 
pairs. The renormalized onsite potential for the regu- 
lar and dimer sites will be denoted as and ed and the 
coupling within the dimers is denoted as 7. It turns out 
that the dimers are connected to each other and rest of 
the lattice by a coupling of unit strength. The renormal- 
ized onsite potential as well as couplings depend upon E 
and are related to the corresponding bare values of the 
parameters by the following equations. 
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Er = {EtEs — 1) 

Ed = {EbEs — 2) 

7 = lEs 



(2) 
(3) 
(4) 



Here E^ = E — ex, where x = s,b, r, d. 

Next we look for a traveling wave solution on the renor- 
malized lattice. We consider Bloch wave solution where 



a propagating wave e 
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at site n with Bloch number 



k undergoes phase shifts as it encounters the dimer de- 
fects. Let flrd be the phase shift between the pure lattice 
sites and the dimers and be the phase shift within 
the dimer. It should be noted that there is no additional 
phase shift between two dimers. Substituting these trav- 
eling wave solutions in the renormalized TBM determines 
the phase shifts. A simple algebra requiring the phase 
shifts to be real ( corresponding to nondecaying wave ) 
determines the resonance condition , 



{E-eb){E 



E = 



(1-7^) 
4cos^(fe7r) 



By eliminating E, we obtain 



(fib - e^) -(1 - 7^)cos(fc7r) 
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(5) 
(6) 

(7) 



The resonance condition we obtain is identical to that 
of Wu et al 1^ for binary alloy. Furthermore, for eb = 
—ts = OL, it reduces to the commutativity condition ob- 
tained by Macia et al ||]. Therefore, we conclude that 
the extended critical states of Macia et al |^ are in fact 
the resonant states. The renormalization analysis |^ dis- 
cussed later on in this paper will provide a confirmation 
of the fact that the resonance condition is globally valid. 
Therefore, for the parameter values at which the reso- 
nance condition is satisfied we have Bloch wave solutions 
on the decimated sites. However, on the original lattice, 
the Bloch waves exist only on the b sites and the wave 
attenuates on s sites and then recovers back on h sites 
whenever resonance condition is satisfied. 

The resonance criterion shows that the off-diagonal dis- 
order, namely 7 different from unity is crucial for obtain- 
ing extended states in the system. It is interesting to 
note that the off-diagonal disorder 7 also determines the 
line of inversion symmetry, E(n,a) = —E{Fn — n^—a) 
for the spectral plot as shown in figure 2. Here Fn is 
the size of the system. This line of inversion symme- 
try coincides with the line where the resonance condition 
is satisfied. The striking feature of this spectral plot is 
the existence of band crossings which form a subset of 
resonant states. It turns out that the states at the band 
crossings corresponds a discrete set of energies with Bloch 
number k — na where n is an integer. Using the expres- 
sion of the transmission coefficient given by Macia et al 
il ^ = (i+c^sin{Ff,7vk)) ' where c is a function of parame- 
ters of the system, we see that these states corresponds 



to T = 1 as the size of the system Fn approaches 00. It 
should be noted that this approach to quasiperiodic limit 
is using Fibonacci sizes only. Therefore, the states at the 
band crossings have a well defined quasiperiodic limit. 
This is in contrast to the rest of the resonant states where 
T oscillates between its minimum value ( corresponding 
to sin{FNTTk) — 1 ) and the maximum value equal to 
unity, as the size of the system changes. 

It turns out that the resonant states at the band cross- 
ings are related to the harmonics of the sine wave when 
expressed in terms of a continuous variable 9 = {na+4>o} 
(where the brackets denote the fractional part and (po is 
an arbitrary constant) except for a discontinuity at a^. 

^^{6) = sin{neTT),0< 9 < 



-sin{Ti9TT), a'^ < 9 < 1 
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(8) 
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In contrast to n-even case, for n-odd, the wave function is 
a double valued function of 9. This analytical expression 
was first observed in our numerical computation of the 
wave function and is shown in figure 3. 

The discrete onsite potential e„ of the Fibonacci sys- 
tem can be viewed as a two-step function of the variable 9 
with a discontinuity at 6* = cr^ where the potential jumps 
from -a to a . The states at the band crossings share the 
discontinuity of the on-site potential. This characteristics 
is true for all resonant states including those that do not 
correspond to the band crossings. Such states are typi- 
cally described by a wave function which is a multivalued 
function of 9 with a discontinuity at cr^. 

The nature of discontinuity in the wave function can 
be understood from our figure 1. As we discussed earlier, 
at resonance, the Bloch wave solutions exist on the deci- 
mated model with r and d type sites after the s sites have 
been eliminated. It is precisely at these s sites that the 
wave first decays in amplitude by 7 and then recovers. 
Therefore, at resonance, although we have Bloch waves 
on the decimated model, the original lattice consists of 
waves that decays but recovers again. Therefore, for in- 
finite systems these waves are always fully propagating 
modes and there is no net attenuation. 

We have described two types of resonant states: the 
states which are related to all harmonics of the sine func- 
tions and have the transmission coefficient equal to unity 
in quasiperiodic limit and the states where the transmis- 
sion coefficient oscillates as we approach quasiperiodic 
limit. We next seek a better characterization of the res- 
onant states using our recently developed RG approach. 
This methodology provides an independent confirmation 
of the existence of extended states in the system. Here 
the extended states are distinguished from the critical 
states by the trivial attractor of the renormalization flow. 
This approach was recently used to demonstrate the ex- 
istence of Bloch wave type phonon modes in the super- 
critical incommensurate Frenkel-Konterova (FK) model. 
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Basic idea underlying the renormalization scheme for 
the quasiperiodic system with golden mean incommen- 
surability is to decimate out all lattice sites except those 
labeled by the Fibonacci numbers F^- Starting from an 
arbitrary initial site n, at the to*'' step, the decimated 
TBM can be written as 

Tpin + Fm+i) = c„(n)V'(n + i^m) + dm{n)'4){n). (10) 

Using the defining property of the Fibonacci numbers, 
Fm+i = F^ + with (Fo = 0,Fi = 1), the follow- 

ing recursion relations are obtained analytically for the 
decimation functions c„i and dm, 

c,„+i(n) = Cm{n + Fm)Cm-l{n + Fm) - (11) 

dm[n) 

drn+l{n) = -drn{n)[dm{n + F.ra) + 

Cm{n -f Fm)dm-l{n + Fm)]c„^(TO). (12) 

These recursion relations can be iterated numerically 
for a large number of decimation steps limited only by 
the precision of the parameters and the energy E and 
the machine precision. For resonant states, since the en- 
ergies are known analytically, the asymptotic behavior 
of the RG flow can be obtained with very high preci- 
sion. Our previous studies have shown that the extended 
eigenf unctions lead to an asymptotic trivial attractors of 
the RG flow while the critical states at the band edges 
are characterized by nontrivial attractors. In the local- 
ized phase, the RG flow decays to zero or becomes un- 
bounded. Therefore, by studying the variation in the 
renormalization attractor as the parameters a and 7 are 
varied, we can determine the nature of eigenstates of the 
system. 

Figure 4 shows the renormalization attractor as a 
varies for a fixed value of 7 for states corresponding to 
the minimum energy (k — 0). There exists a symmet- 
ric 4-cycle which varies continuously as the parameter 
a changes. This implies that the wave function repeats 
itself at every 4*'' Fibonacci site and hence describes a 
translational invariance in Fibonacci space for quasiperi- 
odic system. The variations in the RG attractor implies 
that the scaling properties of the critical wave function 
changes continuously as the parameter a changes. For a 
special value of the parameter which satisfy the resonance 
condition, this nontrivial RG attractor degenerates to a 
trivial cycle ±a~^, determined by the underlying incom- 
mensurability. Therefore, the existence of a crossings in 
the RG attractors signal the resonance transitions. Fur- 
ther studies of resonant states other than those at band 
crossing shows that all such states are described by triv- 
ial attractors of the renormalization fiow. The period of 
these attractors vary as the parameters change but the 
RG flow always settle asymptotically at values which are 
1,0,00 or powers of a. For example, for for E = —5/6, 
the RG flow converges to —1, — ct, —1, 0, 00, 1, 1/ct, 1, 0, ... 



It should be noted that in contrast to the resonant states, 
off-resonant states are always characterized by a non- 
trivial RG attractors. The critical states at the band 
edges are characterized by a period-4 cycle while the rest 
of the critical states are described by strange attractors 
of the RG flow. 

In summary, the RG analysis confirms the fact that 
quasiperiodic Fibonacci systems exhibit Bloch wave type 
solutions. Appropriate decimation scheme traces the ori- 
gin of these states to hidden dimer type correlations in 
the system. Therefore all states that satisfy resonance 
condition are Bloch waves on a decimated lattice. These 
includes the semi transparent states of finite Fibonacci 
lattices discussed previously, (s) For infinite size system, 
these states are fully transparent states with oscillating 
T as there is no net attenuation. Intriguing result of our 
paper is the existence of extended states with Bloch num- 
ber equal to the multiple of the golden mean. In contrast 
to rest of the resonant states, these states which reside at 
the band crossings have transmission coefficient equal to 
unity in the quasiperiodic limit. Therefore, our studies 
provide a clear characterization of extended states in the 
Fibonacci system and shows that the existence of fully 
propagating states is not related to the quasiperiodic long 
range order of the system. 
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FIG. 1. Fibonacci lattice with two types of diagonal cou- 
plings defining a lattice with b and s types of atoms. The 

springs denote the coupling 7 while straight linos denote unit 
coupling. Renormalized lattice after the s type sites are dec- 
imated. The resulting lattice can be viewed as a regular r 
type lattice with d type de wh always appear in pairs. 



FIG. 2. Band spectrum as a function of a for a fixed 7 
equal to 2 in (a) and is equal to unity in (b). The resonant 
states lie on the dashed line which happens to be the line of 
inversion symmetry and passes through the band crossings. In 
the absence of off-diagonal disorder, the resonance and band 
crossings happen at trivial value of a. 

FIG. 3. Figure shows the wave functions at four different 
band crossings as a function of 9: (a) shows the wave function 
for n = 2, 4 which is related to the first and second haxmonics 
of the sine wave, and (b) corresponds to n = 1, 3. Unlike the 
even n case, the wave function is not a single valued func- 
tion of 0. The wave functions show the discontinuity of the 
underlying potential at golden mean. 

FIG. 4. Renorrnalizatioii attractor as a function of a 
showing how the period-4 limit cycle of the decimation func- 
tion Cn changes as a varies for the states corresponding to 
minimum energy(A; = 0) at a fixed value of 7 = 2. The exis- 
tence of resonant states appears as a crossing as the nontrivial 
cycle degenerates to a trivial 2-cycle ±a~^. 
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